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For the phenomenological description of magnetohydrodynamic turbulence competing models 
exist, e.g. Boldyrev [Phys.Rev.Lett. 96, 115002, 2006] and Gogoberidze [Phys.Plas. 14, 022304, 
2007], which predict the same Eulerian inertial-range scaling of the turbulent energy spectrum 
although they employ fundamentally different basic interaction mechanisms. A relation is found 
that links the Lagrangian frequency spectrum with the autocorrelation timescale of the turbulent 
fluctuations, r ac , and the associated cascade timescale, r cas . Thus, the Lagrangian energy spectrum 
can serve to identify weak (r ac <C T cas ) and strong (r ac ~ r cas ) interaction mechanisms providing 
insight into the turbulent energy cascade. The new approach is illustrated by results from direct 
numerical simulations of two- and three-dimensional incompressible MHD turbulence. 



Magnetohydrodynamic (MHD) turbulence has been 
the subject of intense research during the last decades 
since turbulent low-frequency, long-wavelength fluctua- 
tions conveniently described in the incompressible MHD 
framework are present in many astrophysical systems (see 
for example pQ and references therein). As of now, the 
theoretical description of the universal statistical proper- 
ties of turbulent flows relies mainly on phenomenological 
models like Kolmogorov's K41 picture [2] of hydrody- 
namic turbulence. In the incompressible MHD case the 
phenomenological description is complicated by the pres- 
ence of shear Alfven wave modes and their associated 
timescale, the Alfven time ta- Especially for turbulence 
in the presence of a strong mean magnetic field there exist 
competing phenomenological models of the inertial-range 
energy cascade [3HZ]- Although based on physically dif- 
ferent mechanisms the most recent models by Boldyrev 
[6] and Gogoberidze [7] predict identical diagnostic signa- 
tures, e.g., the same inertial-range scaling of the energy 
spectrum, and are thus hardly distinguishable by conven- 
tional measurements in the Eulerian frame of reference. 
Complementing Eulerian diagnostics, the Lagrangian fre- 
quency spectrum gives insight into the timescales associ- 
ated with the turbulent energy cascade and the underly- 
ing nonlinear interactions of turbulent fluctuations. 

This Letter reports a fundamental relation between 
the Lagrangian frequency spectrum and the character- 
istic timescales of turbulence, the autocorrelation time 
r ac and the cascade time T cas . Here, the autocorrelation 
time T ac = T ac (£) characterizes the dominant nonlinear 
interaction process between turbulent fluctuations with 
I being the spatial scale under consideration. On the 
cascade time scale r cas the fluctuations at a fixed spa- 
tial scale loose their coherence and decay into smaller 
turbulent fluctuations. The relation presented in the fol- 
lowing allows to investigate a fundamental aspect of tur- 



bulent dynamics, i.e. to distinguish whether nonlinear 
interactions are weak, r ac <C r cas , or strong r ac ~ r cas . 
The relation is supported by high-Reynolds-number di- 
rect numerical simulations of two- and three-dimensional 
MHD turbulence. In the following, pertinent turbulence 
phenomenologies are summarized focusing on their char- 
acteristic timescales. 

In the K41 picture of Navier-Stokes turbulence the 
auto-correlation timescale r ac is determined dimension- 
ally by the nonlinear turnover time tnl = £/ve where 
vg is a velocity fluctuation (eddy) at scale I. A single 
nonlinear interaction between eddies reduces the coher- 
ence of an involved fluctuation so significantly that it 
ceases to exist at scale £ having generated fluctuations at 
slightly smaller scales. The required time for the associ- 
ated decorrelation in the K41 picture is r cas ~ tnl ~ Tacj 
thus the nonlinear interaction is strong. In incompress- 
ible two-dimensional MHD turbulence the Iroshnikov- 
Kraichnan (IK) phenomenology [3J2] seems to apply, see 
e.g. [5]. There, colliding counter-propagating Alfvenic 
fluctuations lead to a wave-based nonlinear decorrela- 
tion of turbulent structures. The fundamental interac- 
tion timescale is the Alfven-time ta = l/bo where the 
assumption of constant mass density yields the Alfven 
speed as the value of a properly normalized external mag- 
netic guide field bo or, if no such field is present, the 
slowly varying large-scale magnetic fluctuations 6 rms . As 
colliding Alfven wave packets experience only a small de- 
formation in a single nonlinear interaction, many con- 
secutive interactions are required for a cascade step, i.e. 
ta < t nl/ t a ~ r cas since generally r A < r NL , and the 
nonlinear interaction is weak. In incompressible three- 
dimensional MHD turbulence with weak to moderate 
mean magnetic field the Goldreich-Sridhar (GS) [5] phe- 
nomenology has become widely accepted. It enhances 
the IK picture by explicitly taking into account dynam- 
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ical anisotropy with regard to the direction of the local 
magnetic field. In addition, the hypothesis of a criti- 
cal balance of turnover and Alfven time is made, i.e. 
T ac ~ T A ~ T nl ~ T cas resulting in strong nonlinear in- 
teraction. 

However, the inertial-range scaling of the Eulerian en- 
ergy spectrum of MHD turbulence in a strong mean mag- 
netic field does not agree with the GS phenomenology 
A proposed alternative is the dynamic-alignment 
model (DA) [6] which extends the GS model by a scale- 
dependent polarization of the interacting wave packets 
but still leads to strong nonlinear interaction. By a dif- 
ferent approach, an anisotropic variant of the IK-picture 
(AIK) [7] yields weak interaction dynamics. It proposes 
nonlocal decorrelation effects on inertial-range scales by 
large-scale fluctuations. Both phenomenologies lead to 
identical scaling results with regard to Eulerian two-point 
statistics although their underlying physical assumptions 
are fundamentally different. 

The inertial-range scaling of the Eulerian energy spec- 
trum E(k) serves as a standard diagnostic for turbulence 
investigations, but is in fact not unique: the K41 and 
GS models yield E{k) ~ fc~ 5 / 3 . In the GS picture the 
wavenumber k is defined in a direction perpendicular to 
the magnetic field, k — > k±. The DA, AIK and IK mod- 
els, in contrast, all predict E(k) ~ fc~ 3 / 2 with k — > kj_ for 
DA and AIK. This has motivated the development of a 
new diagnostic approach that allows to probe the relation 
between r ac and r cas as presented in the following. 

The Lagrangian two-point two-time velocity correla- 
tion [14] is defined as 

R L = (V(X + r,T + t )-V(X ,t )), (1) 

where V(Xo + r, to + t) is the velocity measured at time 
to + r of a fluid element that was at position Xq + r 
at time to- The Lagrangian variable X — X(Xo,t) 
is the time-dependent position along a fluid particle's 
trajectory. The Lagrangian velocity V(X ,t) is con- 
nected to the Eulerian velocity field by V(X ,t) = 
v(x = X(Xo,t),t). While i?L is generally a correla- 
tion function involving two fluid particles at positions 
XW(X + r,T + t ) and X ( ^(X ,to) it reduces to a 
two-point velocity correlation along a single trajectory if 
r = 0. The corresponding Eulerian correlation function 
is 

i? E = (v(x + r,t + T) -v(x,t)). (2) 

with v(x,t) representing the velocity field at a position 
x and time t while r and t stand for independent trans- 
lations in space and time, respectively. 

For statistically homogeneous and stationary turbu- 
lence these functions depend on r and r only, -Rl.e = 
Rl,e(t,t). The corresponding two-time spectral func- 
tions Q(k, t) are defined as the Fourier transforms of the 



correlation functions 

Rl,e = J d 3 kexp(ik ■ r)Q UE (k,T) . (3) 

The two-time spectral functions are related to the three- 
dimensional energy spectrum E(k) by 

Q hiE (k,r) = E(k)G h , E (r/rlf), (4) 

where r ac ,E are the Lagrangian and Eulerian autocorrela- 
tion time scales, and G^^t/t^, ; ) are the corresponding 
response functions [15] , The general features of mod- 
els of the response function used in the theoretical de- 
scription of turbulence [T5HT7] are that Gl,e is a smooth 
function with Gl,e(0) = 1, Gl,e(x) = for x ^> and 
Jo°° G\ J ^(x)dx = 1. 

The Lagrangian and Eulerian frequency spectra are 
defined as Q3] 

$l,e(w) = ^ J drcos(wT)i? LiE (0,T) . (5) 

Putting Eq. Q in Eq. Q, and plugging the result in 
Eq. ([5]) yields for the frequency spectrum 

$L,E(w) = i Jdr J &kE{k)G^ E {T/ T lf)cos{LOT) (6) 

with E tataX = J dkE(k). The energy spectrum E{k) and 
the associated wavenumber has to be defined in a suitable 
way for the considered geometry e.g. spherical, planar, 
or cylindrical. Under the assumption of self-similarity of 
all involved dependent variables a dimensional approxi- 
mation of this result linking wavenumber and frequency 
energy spectra in a simple way yields 

w*l,e(w) - kE(k) with w-l/r^ E (fc). (7) 

Note that relation ^ is also found in [18] following a 
different and more specific approach. 

The Eulerian correlation time at a fixed position is 
dominated by the sweeping of small-scale fluctuations by 
the largest-scale eddies t^ c ~ (kvo)^ 1 - Consequently, 
the spectral scaling of frequency and wavenumber spec- 
tra should be identical in this case with <i>g S (u;) ~ 
(evo) 2 / 3 oj~ 5 / 3 or <& e k (uj) ~ e 1 / 2 ^^ 3 / 2 where the choice 
depends on the respective Eulerian inertial-range scaling. 
For MHD turbulence, however, the Lagrangian frequency 
spectrum allows to distinguish turbulence phenomenolo- 
gies based on strong nonlinear interaction such as GS and 
DA from pictures based on inherently weak interaction 
like (A)IK. Here, the Lagrangian autocorrelation time is 
assumed to be characteristic for nonlinear interactions 
in the energy cascade. In the MHD case, this holds as 
long as spectral kinetic and magnetic energy are suffi- 
ciently close to equipartition. A detailed investigation of 
the influence of the large wavenumber contribution of the 
Eulerian spectrum on the Lagrangian inertial range |14j 
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TABLE I: Parameters of the numerical simulations. Re: 
Reynolds number, -it lms ; 6 r m S : RMS value of velocity and mag- 
netic field fluctuations; bo: external mean magnetic field; e K , 
e M : kinetic and magnetic energy dissipation rates; v kine- 
matic viscosity; iV co ii oc : numerical resolution. 

shows that the scaling of the Lagrangian inertial range is 
not adulterated by the Eulerian large-scale interval. 

With lo ~ I/t^cO) and E(i) ~ w$ L (w), the classi- 
cal constant-flux ansatz E(£)/r cas ~ e gives ~ 
er cas T^ c . Strong interaction turbulence is characterized 
by T ac ~ r cas resulting in 

$i»~ew- 2 . (8) 

This scaling is thus expected for the GS and the DA 
model and is well-known from Navier-Stokes turbulence 
|2"U] . In contrast, for IK-based phenomenologies 
with a weak interaction mechanism r ac <C r cas ~ T^ L /T ac . 
Using relation r ac ~ uj^ 1 ~ (fcfo ) _1 , and tnl ~ 
(fc 3 £ fc )~ 1/2 Yields r cas ~ b 2 ^ 1 ^ 2 . Thus 

* L (w) ~ £ 1/2 &o^ 3/2 . (9) 

Note that this result holds for IK. In the case of AIK 
the quasi-constant large-scale magnetic field &o has to be 
replaced by w rms . 

Lagrangian frequency spectra are obtained by tracking 
fluid particles in direct numerical pseudospectral sim- 
ulations of MHD turbulence. Details of the numerical 
method can be found in [jJJ [21] ■ In the three-dimensional 
simulations the number of tracers amounts to 3.2 • 10 6 , 
except for the lowest resolution runs where it is lowered 
to 5 • 10 5 . In the two-dimensional simulations 2 • 10 6 
tracers have been tracked. Important parameters and 
characteristics of the simulations are listed in Table [J 
The magnetic Prandtl number Pr m , the ratio of kinematic 
viscosity v to magnetic diffusivity 77, is unity. Significant 
deviations from this value would lead to a degradation 
of observable inertial-range self-similarity due to numer- 
ical resolution constraints and are thus not considered 
here. In the macroscopically isotropic cases, bo = 0, both 
the magnetic and velocity field are forced by independent 
Ornstein-Uhlenbeck processes (221 [23] in the wavenumber 
shell kf = 3 in order to maintain quasistationary turbu- 
lence. In the anisotropic MHD case, b n = 5, large scale 
Alfvenic fluctuations are excited by the stochastic forc- 
ing method which impair the frequency scaling range. 
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FIG. 1: The compensated Lagrangian frequency spectra in 
the macroscopically isotropic three-dimensional (black lines, 
compensated by uj 2 /e) and two-dimensional (grey lines, com- 
pensated by w 3 ^ 2 /(e 1 ^ 2 6 rms )) cases for various Reynolds num- 
bers: 1050 (black dashed), 3150 (black continuous), 2150 
(grey dash-dotted), 6200 (grey dashed) and 18240 (grey con- 
tinuous). The frequency axis is normalized by the Kol- 
mogorov frequency lo v = i^\/e K /v. For clarity the spectra 
for the two-dimensional cases have been shifted by a constant 
factor. 

Therefore in this case the system is forced by freezing 
the lowest wavenumber modes of the velocity and mag- 
netic field of a fully developed turbulent state, see, e.g., 
[llj . In the two-dimensional case turbulence is main- 
tained by keeping the kinetic and magnetic energy in the 
lowest wavenumber shells l<|fc|<3 = fefata con- 
stant value. Based on the estimate (io/^diss) 4 ^ 3 for the 
width of the inertial-range, the Reynolds number is de- 
fined as Re = [27r/(M d iss)] 4/3 where £ diss ~ (i/ 3 /e K ) 1/4 
is the Kolmogorov dissipation length. To obtain the La- 
grangian frequency spectrum first the autocorrelations of 
the velocity and the magnetic field fluctuations along the 
particle trajectories are calculated. The spectrum is then 
computed as the cosine transform of the auto-correlation 
functions 24 

$l,»M = ^JdT((Vi(t + r)Vi(t)) + (bi(t + r)h(t))) cos(u;t) . 

(10) 

In the macroscopically isotropic MHD cases the fre- 
quency spectrum does not depend on the component i 
of the velocity or the magnetic field fluctuations. Here 
the total frequency spectrum $(w) = I2j = i$i(w) ^ s 
shown. In the hydrodynamic case (not shown) the spec- 
trum scales as cu~ 2 in the inertial range in agreement 
with previous experimental [25] and numerical [2U] re- 
sults. In the macroscopically isotropic 3D MHD case 
the Lagrangian frequency spectrum also shows a scal- 
ing with lj~ 2 (see Fig. [l]). This supports the GS phe- 
nomenology with its strong cascade mechanism for this 
configuration. In contrast, the Lagrangian spectra from 
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FIG. 2: The compensated Lagrangian frequency spectrum 
(cf. caption of Fig. [I]) in the MHD case with a strong mean 
magnetic field for Reynolds numbers of about 4410 (cf. table 
ffp. Continuous lines: bo-perpendicular components; dashed 
lines bo-parallel component. 



the two-dimensional simulations which are shown in the 
same figure display an approximate w~ 3 / 2 -scaling indica- 
tive of a weak interaction cascade. This observation cor- 
roborates the IK phenomenology for the two-dimensional 
configuration. In the anisotropic MHD case a depen- 
dence of the scaling on the component of the fluctua- 
tions is observed (see Fig. [2]). For the component along 
bo o;~ 3 / 2 -scaling is observed whereas for the perpen- 
dicular components a scaling exponent a = 1.8 ± 0.1 
(where $(w) ~ uj~ a ) is measured. This suggests that ei- 
ther the energy cascade changes its dynamical character 
from weak interaction (field-parallel) to strong interac- 
tion (field-perpendicular) , or that the energetic structure 
of the flow can not be captured adequately by a simple 
parallel/perpendicular decomposition. 

In summary, a new relation is introduced that re- 
lates the nonlinear auto-correlation time and the cas- 
cade time of turbulence with the Lagrangian frequency 
spectrum. The relation is corroborated by comparing 
high-Reynolds-number direct numerical simulations of 
two- and three-dimensional MHD turbulence with cur- 
rently accepted phenomenological expectations. As the 
Lagrangian frequency spectrum is sensitive to the under- 
lying cascade mechanism it provides additional insight in 
the yet not fully understood case of MHD turbulence in 
a strong mean magnetic field. This is particularly useful 
in cases where the discrimination between different the- 
oretical models of turbulence is hard or even impossible 
to achieve by Eulerian two-point statistics. The ability 
to investigate basic characteristics of turbulent energy 
transfer adds significant value to this approach. 
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